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Abstract. A very well-known machine model in scheduling allows the 
machines to be unrelated, modelling jobs that might have different char- 
acteristics on each machine. Due to its generality, many optimization 
problems of this form are very difficult to tackle and typically APX-hard. 
However, in many applications the number of different types of machines, 
such as processor cores, GPUs, etc. is very limited. In this paper, we ad- 
dress this point and study the assignment of jobs to unrelated machines 
in the case that each machine belongs to one of a fixed number of types 
and the machines of each type are identical. We present polynomial time 
approximation schemes (PTASs) for minimizing the makespan for multi- 
dimensional jobs with a fixed number of dimensions and for minimizing 
the Lp-norm. In particular, our results subsume and generalize the ex- 
isting PTASs for a constant number of unrelated machines and for an 
arbitrary number of identical machines for these problems. We employ a 
number of techniques which go beyond the previously known results, in- 
cluding a new counting argument and a method for making the concept 
of sparse extreme point solutions usable for a convex program. 



1 Introduction 

One of the most general models in machine scheduling is the model of unre- 
lated machines, where the characteristics of each job depend on the machine 
that executes it. As the term "unrelated" suggests, these characteristics might 
be completely different for each machine. This is a very general model which 
causes a significant increase in necessary algorithmic effort and complexity in 
comparison with simpler models. For instance, for makespan minimization there 
are polynomial time approximation schemes known for an arbitrary number of 
identical and uniform machines [617] . but for unrelated machines approximating 
the makespan within a ratio lower than 3/2 is NP-hard. Even more, improving 
upon the best known 2-approximation algorithm |12) has been an important 
open problem in scheduling for more than 20 years. 

In recent years, the design of modern hardware architectures has seen an ad- 
vent of heterogeneous processors in a system, e.g., the cores of a CPU, graphics 
processing units or floating point units. These devices usually have very different 



characteristics since they are especially designed for certain operations. However, 
the total number of different types of processors in a system is usually very lim- 
ited. This motivates the study of unrelated machine scheduling in the setting that 
the given machines (arbitrarily many) are partitioned into a constant number of 
types. This setting subsumes and generalizes two classical scenarios: scheduling 
an arbitrary number of identical parallel machines, and scheduling a fixed num- 
ber of unrelated parallel machines. We study two problems in this setting: first, 
the assignment of multidimensional jobs (in an application, each dimension may 
correspond to a scarce resource such as execution time, memory requirement, 
etc.). The objective is to minimize the makespan across all dimensions. For any 
e > and any fixed number of dimensions, we provide a (1 + e)-approximation 
algorithm for this problem. Our second problem is the minimization of the L p 
norm of the load vector of one-dimensional jobs (instead of the L^ norm). In 
this case we derive a PTAS for any fixed p > 1. 

Related work. Lenstra, Shmoys and Tardos |12| and Shmoys and Tardos |14| pre- 
sented a 2-approximation algorithm for makespan minimization of an arbitrary 
number of unrelated machines. It is known that approximating the same prob- 
lem within a ratio lower than 3/2 is NP-hard [T2]. In fact, the problem remains 
APX-hard for any fixed L p norm [5]. 

Azar and Epstein [2] considered the minimization of the L p norm when 
scheduling unrelated parallel machines. They give a 2-approximation algorithm 
for any fixed L p norm (p > 1) and a PTAS for the case of a fixed number of 
machines. They also give a (D + l)-approximation algorithm for the minimiza- 
tion of the L p norm of the generalized load vector of .D-dimensional jobs, for any 
fixed p > 1. The approximation ratio of the former algorithm was later improved 
to less than 2 by Kumar et al. |10] , 

For the special case of identical machines, a PTAS for makespan minimization 
was given by Hochbaum and Shmoys [6]. A PTAS was later given for general 
L p norms by Alon et al. pQ. For the Z3-dimensional makespan minimization 
problem with identical machines, Chekuri and Khanna j4] provide a 0(log 2 D)- 
approximation when D is arbitrary, and a PTAS when D is fixed. A PTAS for 
makespan minimization on uniform machines has been provided by Hochbaum 
and Shmoys [TJ. 

Our contribution. In this paper we study scheduling problems on unrelated ma- 
chines which are partitioned into a constant number of types, such that two 
machines of the same type are identical. In this setting, we present polyno- 
mial time approximation schemes for the problems of minimizing the makespan 
of £>-dimensional jobs, for constant D, and for minimizing the L p -norm of 1- 
dimensional jobs, for constant p > 1. Both results subsume and generalize the 
known PTASs for an arbitrary number of identical machines |1|4|6| and for a 
constant number of unrelated machines |2I8I12| . Not surprisingly, certain ideas 
in the latter algorithms are useful in our setting as well, for example, geometric 
rounding and enumeration techniques. However, obtaining our results requires 
non-trivial extensions of the known methods since we face obstacles that do not 
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occur in either of the two subsumed cases. In particular, there is no direct way 
to enumerate the assignment of the large jobs, i.e., jobs that are longer than an 
e-fraction of the load of their respective machine in an optimal solution. The 
reason is that a job could be large on machines of one type and small on the 
other and in total there can be a superconstant number of large jobs. 

To remedy this, in our algorithm for makespan minimization we use a linear 
program to assign jobs to slots which we enumerate for large jobs. For rounding 
this (sparse) linear program, we use an iterative rounding approach |9I11| . We 
identify constraints in the LP that can be dropped without affecting the com- 
puted solution too much. To this end, we introduce a new counting argument for 
the number of non-integral variables, which is crucial to the approach and may 
be more generally applicable. This counting argument is one of the novelties of 
our contribution. It was also successfully used in j!3j , which evolved in parallel 
to the research presented here. 

When minimizing the L p -norm, matters are even more complicated. Unlike in 
the PTAS for an arbitrary number of identical machines, we cannot assume that 
all machines will have roughly the same load (on identical machines, this would 
hold after some preprocessing [I]). However, this property is important when 
classifying jobs into "large" and "small". We overcome this obstacle by identifying 
some properties of the optimal solution that we can enumerate in polynomial 
time. Moreover, due to the convexity of the objective function, we cannot use a 
linear program but have to employ a convex program (CP). Unfortunately, we 
cannot assume that the computed CP solution has the same sparseness properties 
of an LP extreme point solution, which is what an iterative rounding strategy 
needs. To address this, we use a trick: we take the computed CP solution and 
formulate a linear program based on it, which we show has to be feasible. We then 
compute a (sparse) extreme point solution of this LP and perform an enhanced 
version of the iterative rounding algorithm for makespan minimization. Hence, 
we show how to make the concept of sparse extreme point solutions usable for a 
convex program, an approach that might be useful in other settings as well. 

2 Problem Definitions 

For both of the problems studied in this article, the task is to assign a set of jobs 
J to a set of unrelated parallel machines M. We let n and m denote the number 
of jobs and machines, respectively. 

Makespan minimization of D- dimensional jobs. The input is represented by a 
positive integer cfj for each job j £ J, each machine i £ M, and each dimension 
d £ {1, D}. The objective is to minimize the makespan, given by 

max max > cf „• , 
deli d\ ieM ^— ' '■> 

where for each machine i the set Jj denotes the jobs assigned to i. In other words, 
each machine i has a load of J2jej- c t j m dimension d, and the objective is to 
minimize the maximum load of all machines across all dimensions. 
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We further assume that the machines are partitioned into a set T = {£\, . . . ,1k} 
of K distinct types, where K is assumed to be some constant. For two machines 
i, i' of the same type, one has cfj = cf, j for each job j and dimension d. 

Lp-norm minimization of one- dimensional jobs. The input is represented by a 
positive integer Cjj for each job j E J and each machine i E M. The objective 
is to minimize 



yjeJi jeJ m 
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We again assume that the machines are partitioned into a set T — {£i,...,£k} 
of K distinct types, for some constant K. For two machines i, i' of the same 
type, one has Cjj = Cj/j for each job j. 



3 Makespan Minimization of Multidimensional Jobs 

We present a polynomial time (1 + e)-approximation algorithm for makespan 
minimization in D dimensions on unrelated machines of at most K types, with 
D and K being constants. 

Let e > and suppose that we are given an instance of our problem. First, we 
establish a binary search framework to estimate the optimal makespan. Hence, 
by suitable scaling, it remains to give an algorithm which either asserts that there 
is no solution with makespan at most 1, or which computes a job assignment 
with makespan at most 1 + e. The general idea of this algorithm is as follows: For 
each machine, we classify job into large and small jobs. With an enumeration 
procedure, we enumerate patterns for the big jobs on the machines. One of the 
enumerated patterns will correspond to an optimal solution. Having guessed the 
correct pattern, the remaining problem is to assign each job either to a slot in 
the pattern (then the job is big on its machine) or to the remaining space of the 
machine (then it is small on its machine). We model this problem as a linear 
program. Given that the LP is feasible (otherwise we know that the enumerated 
pattern was wrong), using an iterative rounding approach, we compute a solution 
with makespan 1 + 0(D ■ e). 

Now we present our algorithm in detail. We call a job j large on machine i if 
there is a dimension d such that cfj > e. Note that in each feasible solution, the 
number of large jobs on each machine is bounded by [D/e\ which is a constant. 
For technical reasons, for each large job j on a machine i we redefine each value 
cfj by setting it to max{cf ■, e 2 /D}. This does not increase the makespan of any 
feasible solution by more than e, as the following proposition shows. 

Proposition 1. Let Ji be a set of jobs on a machine i such that J2je.Ji c tj — ^ 
for each dimension d. Let J s C J t and Jb C Jj denote the small and large jobs 
in Ji, respectively. Then J2je.J a c t,j + ^je.h ma,x i c i.j> e2 /D} < 1 + e for each 
dimension d. 
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Next, we round up each input value cfj to the next greater power of — j^. This 
does not increase the objective by more than a factor 1+e. After this preparation, 
we enumerate the patterns of the big jobs on each machine. Intuitively, a pattern 
for a machine i describes the sizes of the jobs running on i. We call a vector 
q = (q 1 , q D ) a large job type if each q d is a power of 1 + e and e 2 /D < q d < 1. 
Let Q be the set of all large job types. Note that since D and e are constants, \Q\ 
is also bounded by a constant. We call a vector tt G {0,..., [D/e\}Q a pattern 
of large jobs for a machine. For each machine type £, we enumerate how many 
machines follow which pattern. With k being the (constant) number of possible 
patterns for a machine, there are at most G poly(m) combinations for 

the patterns of all machines. Note that since all machines of the same type are 
identical, the actual ordering of the machines of the same type does not matter. 

Assume for the ease of presentation that we correctly guessed the pattern 
which corresponds to an optimal solution. For a machine i with a pattern 7r = 
(7r 9 )g 6 Q, we obtain ir q slots for large jobs of type q and a certain amount of 
remaining capacity iem d (i) in each dimension d. Denote by S the set of all 
slots. It remains to determine an assignment of the jobs to the slots and to the 
remaining capacity on each machine. Of course, a job j can only be assigned to 
a slot s on a machine i if its size on i corresponds to s. Also, a job j can only be 
assigned to the remaining space on i if j is small on i (otherwise we would have 
enumerated a slot for it). We model this assignment problem with the following 
linear program, denoted by Slot-LP: 

(Slot-LP) + Xs 'i = 1 V j G J (1) 

x s ,j < 1 Vs G S (2) 

J2 c tj -a %3 ^ rem d (i) Vi G M, Vd = 1, . . . , D (3) 

x tJ > Vi G M, Mj e J 

x s j > Vs 6 S, Vj G J- 

If Slot-LP is infeasible, then in particular there is no integral solution and the 
enumerated pattern was wrong. Now assume that Slot-LP is feasible. With an 
iterative rounding approach (similar to [9 11]) we round the fractional solution. 
We define LPq to be the Slot-LP. In each iteration t we solve a linear program LP t 
which has the same structure as the Slot-LP, but it will involve only a reduced 
set of machines and jobs. Consider an iteration t. We compute an extreme point 
solution x* to LP t . We say that a job j is fractionally assigned to a machine i 
if x* j G (0,1) and it is fractionally assigned to a slot s if x* • G (0, 1). Using the 
sparsity of extreme point solutions together with a useful counting argument, 
we derive the following lemma. 
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Lemma 1. In x* there is either a machine i which has at most 2D small jobs 
fractionally assigned to it, or a slot s which has at most 2 jobs fractionally 
assigned to it. 

Proof. Denote by I the number of variables which equal 1 and by F the number 
of fractional variables in x*. Let s' be the number of slot constraints which are 
still in LP t and let m! be the number of machines for which there are constraints 
of type ([3]) in LP t . Since x* is an extreme point solution, we have that I + F < 
n + s' + D-m'. Also, it holds that n < I + F/2. We claim that F < 2-s' + 2m'-D. 
Assume on the contrary that F > 2 ■ s' + 2m! ■ D. But then 

n < I + F/2 = {I + F) - (F/2) < (n + s' + D ■ m') - (s' +m'-D)=n 

which is a contradiction. For proving the main claim of the lemma, if each 
machine i had strictly more than 2D jobs fractionally assigned to it and each 
slot had more than 2 jobs fractionally assigned to it, then F > 2 ■ s' + 2m' ■ D 
which is a contradiction. □ 

First, we fix all variables that have integral values in x* . If there is a machine i 
such that in x* there are at most 2D small jobs fractionally assigned to it, we 
remove all constraints of type ([3J for machine i from the LP. This is justified 
since after fixing the integral variables of x* , any solution for the remaining 
variables can violate the constraint ([3]) for i by at most an additional value of 
2D ■ e. The second case of Lemma Q] is that there is a slot s which has at most 
two jobs ji, j2 fractionally assigned to it. Intuitively, we seek an integral solution 
in which either j\ or ji will be assigned to s. We model this by removing j±, j2, 
and s from the instance and by adding a new artificial job jo with the following 
characteristics: 

— jo is allowed to be assigned to any slot where either j\ or ji were allowed to 
be assigned to, 

— for each machine i on which only job j\ but not job ji is small, then cfj 2 :— 
cfj for all dimensions d, 

— similarly, for each machine i on which only job j2 but not job j\ is small, 
then cf j Q :— cf - a for all dimensions d, 

— for each machine i on which both ji and j2 are small, we define cf, - o := 

( x ijJ( x *i,h + X W) ■ c th + ( x ijJ( x tji + X W) ■ c i,h for al1 dimensions d. 

We say that jo subsumes the jobs j\ and j% and that jo disposes the slot s. If 
ji or j 2 were already artificial, we say that jo also subsumes all jobs which were 
subsumed by j\ or j'2, and similarly disposes all slots which were disposed by 
jx ° r J2- Denote by LP t +\ the resulting linear program. Propositions [3J and [4] in 
the appendix show that LP t +\ is feasible (in fact, a solution can be constructed 
from x*). Note that by Lemma [1] we perform at least one of the two procedures 
above and hence, in LPt+\ either the number of machines or the number of jobs 
is strictly less than in LP t . Hence, after at most n + m iterations, we obtain an 
integral vector x which assigns all artificial jobs and all original jobs which are 
not subsumed by an artificial job. 
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Lemma 2. In the vector x, each slot has at most one job assigned to it and 
J2jeJ c tj ' — rem<i W + 2-D • e for each machine i and each dimension d. 
It remains to transform the vector x to a solution for all real (i.e., non-artificial) 
jobs, rather than for some of the real jobs and the artificial jobs. One can show 
that for each artificial job j there is a set Jj of real jobs which were subsumed 
by j and a set Sj of slots which were disposed by j. In particular, we can show 
that for each job j' G Jj there is a feasible assignment of the jobs in Jj \ {j'} to 
the slots Sj. The following lemma proves formally the important properties of 
the artificial jobs. 

Lemma 3. For each artificial job j there is a set of original jobs Jj and a set 
of slots Sj such that 

1. Jj n Jji = and Sj H Sj' = for any two artificial jobs j,j', j ^ j' ; 

2. cf j is a convex combination of {cfj,\j' G Jj} for each machine i and each 
dimension d, 

3. no slot in Sj is used by x, 

4- no job in Jj is assigned by x, and 

5. for each j' G Jj there is a feasible assignment of the jobs in Jj \ {j 1 } to the 
set of slots Sj. 

Using the above properties of the artificial jobs we transform x to an integral 
solution for all real jobs. Let j be an artificial job assigned to some slot s £ Sj. 
Then there is a job j' G Jj that can be assigned to s and we assign all jobs in 
Jj \ {j'} to the slots Sj. For replacing the artificial jobs which are not assigned 
to a slot but to the remaining space on a machine we consider all those jobs on 
a machine at the same time. 

Lemma 4. Let i be a machine and let RJi and AJi denote all real and artificial 
jobs, respectively, which were assigned to the remaining space on i. There is an 
integral assignment x' of the jobs in Uj^AJiJj to i and the slots in Uj^AJtSj 
such that each slot gets at most one job assigned to it and 

E c ii • *u + E c t ■ <j • • 

Proof. A solution to the following linear program can be extracted from x: 

(Art-LP), x i,i ^ 1 V -?" e AJ * 

E r t ••'•<•, • E cl r x t , 3 <Tcm. d {i)+2D-e Md G {1, . . . , D] 

Xij > Vj G Uj> e AJiJj>- 

In an extreme point solution to (Art-LP) i there can be at most \AJi \ + D non- 
zero entries. Hence, by a standard counting argument one can show that there 
are at most D sets Jji from which a job is fractionally assigned. Rounding up 
all fractional values and assigning the other jobs in the sets Jj to the respective 
slots in Sj yields a solution with the claimed properties. □ 
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Applying the procedure of Lemma U to each machine i G M yields an integral 
assignment of jobs with a makespan of 1 + 3De in each dimension. Together with 
the binary search framework, this yields our main theorem of this section. 

Theorem 1. Let D.K G N be constants. For any e > there is a (1 + e)- 
approximation algorithm for makespan minimization in D dimensions on unre- 
lated machines of at most K types. 

4 Lp-norm Minimization 

In this section, we present a (1 + e)-approximation algorithm for assigning jobs 
on machines with K types to minimize the L p -norm of the loads of the machines, 
for 1 < p < oo. Note that since we work with the L p -norm, we assume the jobs 
to be 1-dimensional. 

Our strategy is the following: First, we reduce the complexity of the problem 
by enumerating certain structural properties of the optimal solution. Those will 
include the patterns for the big jobs (like in Section [3]) but also certain other 
information which are only important when working with L p -norms. For the 
remaining problem we will formulate and solve a convex programming relaxation. 
Unlike for linear programs, we cannot assume the obtained CP solution to be 
sparse (like e.g., an extreme point solution). However, from the obtained CP- 
solution we will derive a feasible solution to a linear program. From this LP we 
obtain an extreme point solution which we can round by enhancing the iterative 
rounding scheme presented in Section [3] Hence, in our algorithm we make the 
concept of (sparse) extreme point solutions of LPs usable for a convex program. 

Let e > and p > 1 be constants. First, instead of minimizing ||<7|| p , where g 
denotes the vector given by the loads of the machines, we minimize \\g\\ p . Note 
that a (1 + e) p -approximation algorithm for the latter translates to a (1 + e)- 
approximation for the former. Suppose that we are given an instance of our 
problem. We start by enumerating certain properties of the optimal solution. 
In an optimal solution, there might be some machines which execute only one 
job. Intuitively, these jobs are quite large. We call those machines huge. For each 
type t, we enumerate the number of huge machines, denoted by hi. Note that 
since all machines of a type are identical, it does not matter which exact machines 
are huge. Hence, the total number of combinations we need to enumerate is 
bounded by m . For each type £, denote by M^i the huge and by M Sj ^ the 
non-huge machines of this type. For each type t we enumerate the /(p, e) largest 
huge jobs which are processed on a machine of type £, where f(p, e) is a value 
which we obtain from the following proposition. We call them the very huge jobs 
and the corresponding machines the very huge machines. 

Proposition 2. For each p > 1 and each e > there is a number f(p,e) such 
that E ?eG 9 P + (2 • min{ 3 G < (1 + e) p • E g eG 9 P for any set G of positive 

reals with \G\ > f(p,e). 

Knowing the very huge jobs for each type will imply later that we can afford 
making certain mistakes when assigning the remaining huge jobs. Note that there 
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are at most n K '^ p,< ^ possibilities to enumerate. Also, for each type £ we guess 
the longest job which is scheduled on a machine of type £ and which is not huge. 
Denote by c max ^ its length. There are at most n K possibilities for this. 

Lemma 5. Consider an optimal solution and let i be any machine of type £ 
that is not huge. Then its load is at least c maXi £. Moreover, the loads of any two 
non-huge machines of type £ differ by at most c max /. 

Proof. The claim follows from the strict convexity of the objective function and 
an exchange argument, see Appendix lAl for details. □ 

For any machine i of type £ we call a job j huge on machine i if Cjj > c maX) ^. 
For each type £ denote by Hg the jobs that are huge on machines of type £ and 
which are not longer than the f(p, e) very huge jobs for type £ which we guessed 
above. Finally, for each type £ we guess a value ag 6 {1, n) such that in the 
optimal solution the load of each machine of type £ is at least ag ■ c maXj ^ and at 
most (ag + 2) • c max ,£. Due to Lemma [5] such a value ag must exist. Note that 
there are at most n K possibilities for ag. Next, we enumerate the patterns of the 
big jobs on the non-huge machines of each type. To this end, we define a job j to 
be large on a machine i of type £ if cy > e ■ ag ■ c max ^ and small otherwise. Like 
in Section [3] we enumerate over the (polynomial number of) patterns for each 
type. From now on, assume that we know the correct values for all enumerated 
quantities. 

With this preparation, we formulate the remaining problem as a convex pro- 
gram, which we denote as Slot-CP. Like in Section [3j denote by S the set of 
slots for the big jobs. If a job j fits into a slot s then we introduce a variable 
x s j which indicates whether j is assigned to s. For each combination of a job j 
and a machine i such that j is small on i, we introduce a variable Xij. Finally, 
if a job j is huge on machines of type I (i.e., cy > c max ^) then we introduce 
a variable xgj indicating whether j is assigned to one of the huge machines of 
type £. For each machine i, let Bi denote the total length of the large jobs on i. 
Let M s := UgM St g denote the set of all machines which are not huge. For each 
type £, denote by M v h/ the very huge machines. For any very huge machine 
i G M v h t g we define a constant t* denoting its load (due to its guessed job) and 
for any machine i £ M s we introduce a variable t, which models its load. We 
solve the following convex program to an additive error of e. This can be done 
in polynomial time since we have a separation oracle and the objective is convex 
and differentiable [5]. 
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(Slot-CP) min^(i i + ^) p + ^ E K) ? + EE^'M P 

E ^ J + E + E Xi 'i = 1 V J e J ( 4 ) 

E ^.j ^ ^ V/eT (5) 

jeJ 



jeJ 



x i,j — ^2 


Vi G M s 


Cmax ^ 


W g T,V« g M S; . 


Xij > 


V« G M s ,Vj G J 


x s j > 


V.s G 5,Vj G J 


X£j > 


W G T,Vj G J 


u >0 


Vi G M s . 



Since Slot-CP is a relaxation of the original problem, its optimal value yields 
a lower bound on the optimum. Denote by t* the values obtained for the t,- 
variables from the optimal solution of (Slot-CP). In order to be able to use the 
concept of extreme point solutions, we derive the following linear program where 
all the i*'s are constants. 



(Slot-LP) min ^(cejY ■ x t j (7) 

E Xl >i + H Xs '3 + H Xl >i = 1 v i e J ( 8 ) 

x e,j < h e W G T (9) 

jeHe 

J2 x s,j< 1 VsgS* (10) 
jeJ 

E ••''- ./ ' Vi G M s (11) 

jeJ 





> 





Vi 


G M s ,Vj G J 


X S J 


> 





Vs 


G S.Vj G J 


X t ,j 


> 





W 


G T,Vj G J. 



We devise an iterative rounding algorithm which computes an integral solution 
whose overall value is only by a (1 + 3e) p -factor bigger than than the optimal 
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value of (Slot-CP). Like in Section [3] we work with linear programs LP t where 
LPq is the Slot-LP and each LP t is obtained by taking LP t -\ and fixing some 
variables and removing some constraints. In each iteration t, we compute an 
extreme point solution x* of LP t . 

Lemma 6. Let x* be an extreme point solution for the linear program LP t for 
some iteration t. Then either 

1. there is a machine i with at most two small jobs j such that x* ■ S (0, 1), or 

2. there is a slot s with at most two jobs j such that with x* ■ £ (0, 1), or 

3. there is a type £ with at most two jobs j G Hi such that x\ ^ £ (0, 1). 

Proof. We can follow a similar argumentation as in Lemma [TJ via the total 
number of non-integral variables. □ 

First, we first fix all variables which are integral. If either case 1 or case 2 of 
Lemma IS] applies we do the same operation as in Section [31 i.e., drop a constraint 
of type ([lip or replace two jobs by an artificial job and drop a constraint of 
type (|10p . If case 3 applies, i.e., if there is a type I with at most two jobs j S Hi 
such that £ (0, 1), then we define a schedule for the huge machines of type £ 
by assigning each integrally assigned huge job on a single (huge) machine and 
assign the two fractionally assigned huge jobs together on one of the machines 
of type £ which we defined to be huge. We call the latter machine the improper 
machine of type £. Then we remove the constraint © for type I. We will show 
later that the cost of the improper machine is very small in comparison with the 
cost of the very huge machines of the respective type. After the last iteration, 
we replace the introduced artificial jobs by the original jobs that they subsumed 
(see Lemma 21). 

For each machine i, let gi denote its load in the computed integral solution. 
With a similar reasoning as in Section [3J we can show the following lemma. 

Lemma 7. For each small machine i G M s i it holds that gi < t*+3e-a£-c max e < 
(l + 3e)t?. 

Apart from the improper machines, the cost of the huge machines does not 
increase due to our rounding scheme, as the next lemma shows. 

Lemma 8. Let x* be an optimal solution to the Slot-LP and let £ be a type. The 
cost of its huge machines is bounded by 



E s!< E K) P + ( 2 ' min {^ e M ^AT + E x h ■ M p - 

Finally, we show that the cost of the improper machines is small in comparison 
to the cost of the very huge machines. This follows from Proposition [5] and 
Lemma [5] 

Lemma 9. Let £ be a type. Then E iS m hf (**) P + ( 2 ' min{tf|i £ M vh . e }) p < 



ii 



Using Lemmas [7J HI and [9] one can show that the total cost of the final solution 
is at most by a factor (1 + 3e) p larger than the total cost of the optimal solution 
of Slot-CP (see Lemma [TT1 in the Appendix). This yields our main theorem. 

Theorem 2. Let K e N and p > 1 be fixed. For any e > there is a (1 + e)- 
approximation algorithm for assigning jobs to unrelated machines of K types to 
minimize the L p -norm of the load vector. 
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A Appendix 

Proof (Proposition^. By the pigeonhole principle there can be at most D/e 
large jobs in Ji. Thus, 

max{< 3 , e 2 /D} £ cf tj < (D/e) ■ e 2 /D = e 

and thus c i,j + Sje/b ma - x i c i,j, e2 /D} < 1 + e for each d. □ 

Proposition 3. After removing a machine i and the variable Xij for each job 
j , the new linear program LPt+i is feasible. 

Proof. Let ji be a removed job. We decrease x*, - t and x* ^ to zero for any 
remaining machine i' or slot s. This cannot violate any residual constraint. □ 

Proposition 4. Let s be a slot and suppose that in LP t for any job j with j ^ j\ 
and j ^ j2 it holds that x s j £ {0, 1}. After replacing j\ and j% by the artificial 
job jo and removing s, the new linear program LPt+i is feasible. 

Proof. If x* is feasible for LP t , consider the solution x defined as: 







j + jo 


* 

X i,3i 


1 * 

+ X i,h 


j =jo, 






j + 30 


X 8,jl 


4- T* ■ 


j = jo 



It is not hard to check that due to the construction, constraints ^ and of 
LP t +i will be satisfied by x. For example, for constraint (|3|) and the case that 
ji and j2 are both small on machine i, 

5Z c i,i^ij — C i,j X i,3 + C i,j'o( X ji +^2) 

Ed * , d *,c£ * 
C i,j X i,j ~*~ C i,h X 3l + C hh X 32 

< rem d (i). 

Constraint ([T]) may be violated by jo in the sense that Xi j0 + ^ s i s ,j ma y 
be larger than one; however, in that case it suffices to scale uniformly down the 
values (xi JO )igM, i x s,j )seS to obtain a feasible solution. □ 

Proof ( of Lemma\3\) . To obtain the sets Jj and Sj we trace back, in the execution 
of the algorithm, which jobs have been subsumed by j, which jobs have been 
subsumed by those jobs, and so on. Associate to the execution of the rounding 
algorithm a bipartite graph G whose nodes are the jobs and slots that have 
been removed by the algorithm. Whenever two jobs j\, ji competing for slot 
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Fig. 1. Example of the construction in Lemma [3] Gray circles represent original jobs, 
white circles represent artificial jobs, and squares represent slots. 



s are subsumed by an artificial job jo, insert the arcs (jx,s), (j2,s), (s, jo) in 
G. The resulting graph will be a directed forest in which every original job is a 
leaf, every slot has indegree two, and every artificial job has indegree one (see 
Figure Q] for an example). We define Jj to be the leaves of the tree of G having j 
as the root, while Sj is the set of slot nodes of the same tree. Properties (l)-(3) 
are now clear from the definition of the algorithm. Properties (4) and (5) follow 
from the fact that the slots in Sj and jobs in Jj have all been removed from the 
linear program. Property (6) (reminiscent of the fact that, in a tree in which 
every internal node has two children, the number of internal nodes is equal to 
the number of leaves, minus one) is proved by induction on the structure of the 
tree: consider a slot s having two jobs ji, j<x G Jj as children, find inductively a 
feasible assignment for the smaller tree where ji, ji and s have been removed 
and the artificial job that is the parent of s has been replaced with whichever of 
ji; fits jo's parent, and assign the other job to s. □ 

Proof (of Proposition^. Take f(p,e) = (1 + 2f)/(l + e) p . Then 
gp^m^) < (||pIC + 2pM)" P 



1 + 2 p\ 



Ml, 



\G\ 

<(l + e)|| 5 || p . 

Hence, raising both sides to the power of p, we obtain the claim. □ 

Proof (of Lemma\^) . The claim will follow from the convexity of the objective 
function and an exchange argument. Consider the first part of the claim; the 
second part is proved similarly. Let i 1 be any machine of type I where a job of 
length c max .£ is assigned, and let oq, at be the loads of machines i, i' respectively; 
by construction a\ > c max .£- Assume by contradiction that ao < c max ^. Consider 
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the assignment where we exchange all the jobs on i with the job on i' of length 
c m ax,£- Since the machines are of the same type, the processing times of the jobs 
are unaffected. Call a' , a' x the new loads of i,i' . Observe that ao < a' < a\ and 
do < a'i < <Xx] equivalently, there exist /i, rj £ (0, 1) such that a' = /iao + (l— 
a[ = r\a\ + (1 — ?7)ao. However, since a' + a[ = ao + a\ and ao < a%, one has 
fi = r\. Then (a^a^) ~ /i ■ (a ,ai) + (1 — fi) ■ (a l5 ao). Consequently, if a and 
a' are the load vectors of all machines before and after the exchange, and a" is 
the load vector obtained from a by exchanging the i-th and i'-th component, 
a' = fia + (1 — fi)a". Finally, from the strict convexity of the L p norm (for 
1 < p < oo), 

\W\\ P = ha + (1 - ria% < fx \\a\\ p + (1 - \\a% = \\a\\ p , 

which contradicts the assumed optimality of the solution. □ 

Lemma 10. Let T' be the set of all types £ for which there is still the constraint 
^2j eHe xij < hg in LP q and LP q+ \. Then 

E E x h ■ m p E E »ff • m p 

eeT' j£H e e £ T' j£H e 

where x q and x q+1 denote optimal solutions for LP q and LP q +±, respectively. 

Proof. According to the proof of Proposition [31 when in our iterative rounding 
routine we drop a constraint of type the optimal solution for LP q yields 
a feasible solution for LP q+ i. In particular, the value of the optimal solution 
does not increase. By the proof of Proposition 2] the same holds when we replace 
two jobs by an artificial job. Also, when we assign jobs to the huge machines of 
some type £ then the remaining variable assignment stays feasible. Hence, the 
optimal objective value for LP q +i is upper bounded by the optimal objective 
value for LP q . This implies the claim of the lemma. □ 

Proof (of Lemma\S\). For the very huge jobs the contribution is clear. For the 
huge jobs excluding the very huge jobs using Lemma[lO]and the definition of our 
operation for the removing the constraints from Inequality @ we get a bound 
of (2-min{i?|ie AW}) p + E 3 4r M"- D 

Lemma 11. For the loads gi of the computed solution it holds that 

UWp < (1 + 3e) p OPT CP < (1 + 3e) p OPT p 

where OPTcp denotes the value of an optimal solution of Slot- CP and OPT p 
denotes the value of an optimal integral solution of Slot- CP. 

Proof. Recall that the vector t* indicates the loads of the very huge machines 
due to the guessed very huge jobs and the loads of the small machines due to 
the allocation of the small jobs by Slot-CP. Then 
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IK<E E (t*+3e-^- Cmax ,,r+£(i+^ e (si) ? +EE^'fe) p 

e£Ti£M 3: e IdT i£M vh , e t.ET j 

<E E ^+ 3e - < *) p + E( 1 + £ ) p E fe) p +EE4rfe) p 
<(i+3 e ) p (E E ***+E E (*iT + EE^-(^) p 

< (1 + 3e) p OPT C p 

< (1 + 3e) p OPT p . 
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